function cost = C_s(s_dot)

global gamma_s
global delta
global alpha
global c_s_1
global c_s_2
global s

state_dependent_term = max(gamma_s-s,0);

cost = c_s_1*max(s_dot+delta,0) + c_s_2*max(s_dot+delta,0).^alpha+state_dependent_term.*max(s_dot+delta,0);

end

